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Abstract 

Dynamical mass generation in a three-dimensional version of finite-temperature 
QED is studied with the help of Schwinger-Dyson equations in the real-time 
formalism. We go beyond the bare-vertex approximation and include wave- 
function renormalization effects. This introduces a system of two integral equa- 
tions which are solved numerically. In order to increase the reliability of the 
results, fermion and photon self-energies varying independently with energy 
and momentum are used. The method applied enables a detailed study of the 
behaviour of the theory with increasing temperature and number of fermion 
flavours. 

PACS: ll.10.Wx, 11.15.Tk, 12.20.Ds, 12.38.Lg 



* e- mail : georg® phy sik . t u-muenchen . de 







1 INTRODUCTION 



Three-dimensional QED can exhibit interesting non-perturbative phenomena like 
dynamical mass generation JT|-[[|. A finite-temperature version of the theory, re- 
ferred to sometimes as T3 — QED, has been suggested to be relevant to high- 
temperature superconductivity Q . The critical behaviour of the theory with varying 
temperature T and number of fermion flavours Nf was recently studied numerically 
in the Schwinger-Dyson formalism ||. In that paper, a bare- vertex approxima- 
tion was used in the gap equation for the fermion self-energy £, and wave- function 
renormalization was neglected. These approximations correspond to a particular 
truncation of the Schwinger-Dyson hierarchy and to the simplest way of making 
this scheme useful for computations. Since this could influence the final results, we 
return to the problem in order to study the impact of these approximations. Their 
potential importance has been stressed in ||, especially in connection with the be- 
haviour of the theory when the number of fermion flavours is varied. This has led to 
several interesting studies at zero-temperature trying to include these effects JjJ-||, 
in order to solve the controversy on whether there is a critical number of fermion 
flavours beyond which there is no mass generation. At finite temperature, similar 
studies have so far used several approximations that we will relax in this study. 

The first goal of this paper is to provide an estimate for the ratio r = 
2M(0, 0)//cbT c that is free from severe approximations used frequently in the 
past. Here M(0, 0) is the fermion mass function at zero momentum and energy and 
T c is the critical temperature beyond which there is no mass generation. Apart from 
being theoretically very interesting, this quantity is of direct physical relevance since 
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it can be compared to corresponding values measured in certain high-temperature 
superconductors. The second goal is to provide a reliable phase diagram of the 
theory with respect to T and Nf. In order to achieve these goals, the bare- vertex 
approximation is here relaxed and wave-function renormalization is taken into ac- 
count. 

This complicates the study considerably and constitutes a non-trivial step 
forward, since instead of having only one gap equation, a system of two integral 
equations has to be solved. In particular, the numerical algorithm becomes more 
complicated than before and each iteration requires more computing time. As in 
H however, we continue to work with fermion and photon self-energies varying 
independently with energy and momentum, something which by itself is a very 



important improvement of similar investigations in the past [1C]. We continue to 
neglect the imaginary parts of the photon polarization functions and of the fermion 
self-energy for simplicity, as other studies have done so far fj"0[| - [|l2"| . A more detailed 
treatment would of course involve the full propagators , but this is left to future 
studies due to the limited numerical accuracy of the present method. 

2 TAKING WAVE-FUNCTION RENORMALIZATION 
INTO ACCOUNT 

The study of the fermion self-energy and mass function in a non-perturbative context 
is central to our study, and the Schwinger-Dyson formalism will provide the basis 
of this investigation. The Schwinger-Dyson equation for the fermion self-energy is 
given by 



S-\p) = S^{p) - e 2 / ^-fS{q)^{k)VM (1) 



where q = p — k, e is the dimensionful gauge coupling of the theory which we will 
take to be constant throughout this paper, A^ u is the photon propagator with //, v = 
0, 1, 2, T u is the full photon- fermion vertex, 7^ is a four-dimensional representation 
of the 7-matrices, Sq is the bare fermion propagator, and the finite-temperature 
fermion propagator in the real-time formalism is given by 

1 2it5{(l+A(p)) 2 p 2 + Y 2 {p)) 



5(p) = ((l + A(p))ji + E(p))x 



{1 + A(p)) 2 p 2 + T, 2 (p) e^M + l 

(2) 



where (3 = 1/ksT, A(p) is the wave-function renormalization function, 5 is the 
usual Dirac function and we have made a rotation to Euclidean space. Note that 
we avoid the matrix form that the propagator has in this formalism, since the 
Schwinger-Dyson equation that we have written down involves only a one-loop dia- 
gram directly, so complications due to the field-doubling problem do not arise |14fl . 
Moreover, due to the broken Lorenz invariance at finite temperature, the wave func- 
tion renormalization could in principle affect differently the po and \p\ propagator 
components, i.e. we should replace (1+A(p))^> by ({l+A{p))^+a)po+(l+B(p))^ l pi 
with i = 1,2. For simplicity we will restrict ourselves to situations where a = 0, 
which correspond to a zero chemical potential, and we will work in the approxima- 
tion where A(p) = B{p) also for non-zero temperatures as done in similar studies 
| |12| | for simplicity, even though there is no a priori justification for this. We will 
return to this matter later. 

For the vertex T u (p,q) = T a{p, c L)1 v we use two different ansatze, where we 
take T^p, q) to be equal to l + A(q) or to the symmetrized one 1 + (A(p) + A(q)) /2, 
and A(q) is the same wave-function renormalization function appearing in the 
fermion propagator. Even though these vertices do not satisfy a priori the Ward- 



Takahashi identities, they are expected to incorporate the basic qualitative features 
of a non-perturbative vertex at zero temperature when used in a Schwinger-Dyson 



context [15|. The first choice has been used in a finite-temperature case [12], sup 



ported by the qualitative agreement of the results of this ansatz with the ones ob- 



tained by more elaborate treatments ||. The second vertex is discussed in Ref , fig] 
and we will refer to it as the symmetric vertex. Even though these choices corre- 
spond to specific truncations of the Schwinger-Dyson hierarchy, they expected to 
give qualitatively reliable results. 

Moreover, the photon propagator in the Landau gauge is given by [|j 

A ^ {k) = k 2 + n L (k) + k 2 + n T (k) (3) 



where 



k 2 

Qiiu = (Sfj.0 - k^,ko/k 2 )^-{5 U Q - k u k /k 2 ) 

k 

Pvv = 5^(6^ -kikjjk 2 )b v j (4) 



with i,j = 1,2, and where we neglect its temperature-dependent delta-function 
part since it is expected to give a vanishingly small contribution ||, [11], [17|. The 



longitudinal and transverse photon polarization functions Hi and Ht are given in 



[11 1, where they are calculated in a massless-fermion approximation. Note that in 
principle one should also have an integral equation for the polarization functions 
coupled to M and A, but this would complicate this study even numerically too 
much, and we therefore continue to work with this approximation. 

Identifying the parts of this equation with the same spinor structure, we 
reduce the problem to that of a system of two three-dimensional integral equations 



involving two functions varying independently with po and \p\. The equations take 
the following form: 

a f dk \k\d\k\d9 M{q , \q\)R(p , \p\,Qo, \q\) 



I -h\ a f dk k a 



N f (l + A( P0 , \p\)) J (2vr)3 9 2 + M 2 (g0; 



V 1 



a f \k\d\k\dO M(E,\q\)R(p ,\p\,E,\q\) v 



f\k\d\h 

J (2vr 



N f (l + A(p , \p\)) J (2vr) 2 2E{eP E + 1) 



E E 



V-i p=^,t (Po " e^) 2 + ^ + n P (p - e£, \k\) 
a f dk \k\d\k\d9 R(p ,\p\,q ,\q\) 



N f plJ (2vr)3 q 2 +Ml(q ,\q\) ' 
Q(p ,p,k ,k) P(p ,p,k ,k) 



x 



v £; 2 + n L (fco,|fc|) fc 2 + n T (^,|fc|), 

a /■ |jfe|d|fc|d0.R(po, |j5l,£,|<?1) 



iV> 2 7 (27r) 2 2£(e^ + 1) 



x 



y- ( Q(Po,P,Po ~ eE,k) ^ 

=T"-i V (PO " e^) 2 + & + n L (po - eE, k) 



+ P(po,p,po - eE,k) \ 

(p - eEf + P + U T {p - eE, k)J ' 

where a = e 2 iV/, R(po,\p\,qo,\q\) = ^^^"^ , and it is more convenient to 

work with the mass function M(po,\p\) = T,(po,\p\)/(l + A(po,\p\)). The function 

R is obviously equal to the identity for the first vertex choice. We also sum over 
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the photon polarizations P = L,T and over the two roots of the delta function 
by introducing e = 1,-1. The quantity E is approximated by the relation E 2 ~ 
\q\ 2 + M 2 (0,0) pi, where use of the delta-function property 6 (ax) = S(x)/\a\ has 
been made. 

Furthermore, the functions Q and P are given by 

Q(Po,P, k , k) = 2 \p Q p— Jpylo p - J ~PQ 

P(p ,p,k ,k) = 2 Uq-^^\- Pq . (6) 

One can easily check that Q + P = —2(pk)(kq)/k 2 , which would reproduce the 



result of [12| if one takes Iii(k) = Iix(k) = n(&) and switches to imaginary-time 
formalism (note that our q has the opposite sign than the one in |l2| ). 

The equations for M(po, \p\) and A(po, \p\) have to be solved self-consistently, 
and the equation for the mass function actually always accepts, apart from the 
solutions we will seek, the trivial solution as well. An analytical study of such a 
system would not be possible without severe approximations, so we proceed to the 
numerical solution of the equations given above. 

3 SOLVING THE SYSTEM OF EQUATIONS 

3.1 The algorithm 

In order to attack the problem numerically, we have to discretize our momentum 
space. The technique will in principle be similar to the one in ||, but the algorithm 
will be substantially complicated due to the fact that we now have to solve a system 
of two four-dimensional integral equations involving functions of two variables. The 
gap equations are both ultraviolet (UV) and infrared (IR) finite, with a providing 
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an effective UV cut-off [|| and the fermion mass function providing a physical IR 
cut-off. Therefore, we use UV and IR cut-offs Am and Auy only for numerical 
reasons. However, we are always careful to keep the temperature T and the mass 
function M within the range of Am and Ajjy. It is a well-known problem however 
that the existence of this IR cut-off does not allow us to draw firm conclusions 
on the criticality of the theory, since a mass function dropping below this cut- 
off could be small but still not exactly zero. We will return to this issue later. 
The complexity of the algorithm does not allow us to consider Am/Auv ratios 
smaller than 10~ 4 . Since the momentum space spans several orders of magnitude, 
we use logarithmic variables and discretize the squares of external and loop momenta 
according to log 10 (A| iJ ) + ^ log 10 (Ajyy/A^), and the integration angle 6 according 
to 2iri/n, where i = 1, n. We have therefore a five-dimensional lattice, with three 
dimensions coming from the integrals and two from the external momenta. 

The results presented in this paper are for a lattice having 16 points in each 
of the five dimensions. We checked the stability of our results for lattices of other 
sizes as well, although computing-time limitations did not allow us to consider 
much larger lattices. We then use a relaxation method described in Ref.||] to solve 
our system of equations self-consistently by iterations. We do not exceed a 10% 
numerical accuracy of the results in the worst cases, which we judge as satisfactory 
due to the complexity of the algorithm. We trace back this numerical uncertainty 
to the minimum achievable root-mean-square deviation of the input and output 
configurations of the algorithm. This deviation is much smaller if one neglects 
the high-momentum regions where the solutions are vanishingly small and where 
numerical errors appear easier. Note however that the mean relative difference of the 
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input and output configurations is monotonically decreasing with each iteration, like 
the deviation, but reaches values as low as 10 -4 , which we consider as an indication 
of an overall convergence. On the contrary, in non-converging cases this difference 
is increasing and the configurations tend rapidly to zero, apart from the fact that 
the deviation is oscillating randomly taking sometimes unusually large values. We 
therefore have a clear-cut convergence criterium to enable us to study the critical 
behaviour of the theory. 

3.2 Behaviour of the theory with temperature 

The solution at T = and Nf = 2 for the functions S(po, \p\) and —A(po, \p\) for the 
vertex with Ta(p, q) = 1 + A{q) is given in Figs. 1 and 2 respectively. The general 
form of these functions does not change with increasing temperature, varying of Nf 
or by using the symmetric vertex instead. The function S(po, \p\) falls as expected 
with increasing momentum, and is of the same form as the mass function M(pq, \p\). 
The function A(po, \p\) is always in the range between -1 and as required [Iq ], it is 
tending to zero for increasing momenta, and it is of the same form and magnitude 
as the approximate form used in jl^]. In view of the approximation A(p) = B(p) for 
the wave-function renormalization stated at the beginning, the roughly symmetric 
behaviour of the two functions A and X with respect to po and \p\ argues in support 
for the self-consistency of the calculation. Moreover, the low-momentum values 
of the wave-function renormalization, which are on the order of ^4(0, 0) ~ —0.3, 
indicate the importance of including its effects in a complete calculation in this 
gauge. 

For a given number of fermion flavours Nf, when the temperature exceeds 
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Figure 1: The fermion self-energy E(po; \p\) a t zero temperature and for Nj = 2, 
Kjjy/a = 0.1, as a function of energy and momentum in logarithmic scale. All 
quantities are scaled by a. 
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Figure 2: The opposite of the wave- function renormalization — A(po, \p\) at zero 
temperature and for Nf = 2, Kuy/a = 0.1. 
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some critical value there is no solution for the fermion mass function but the trivial 
one. We infer that not only by the fact that the relaxation algorithm does not 
converge, but also because the output configurations tend rapidly to zero and after 
a number of iterations they fall below the IR cut-off. One should note that the 
temperature at which this happens could be somewhat smaller than the critical 
one, since the zero-momentum mass function at slightly smaller temperatures is not 
very much smaller than the one at zero temperature. The abrupt fall of the solution 
with temperature thus could possibly be an indication that the initial smooth input 
configuration is unable to converge to the correct solution, which is rapidly varying 
but still not zero. Since the dimensionless ratio r of (twice) the zero-momentum 
and zero-temperature fermion mass function divided by this critical temperature is 
important for the study of solid-state systems, we list in Tables 1 and 2 its values 
for various choices of Nf and Kuy/a, for our two different choice of vertices. Note 
that the case Nf = 2 is the one relevant to high-temperature superconductivity. 

The fermion mass function falls with increasing momenta and the results 
obtained should not depend on the choice of the UV cut-off. Even though the 
limited number of lattice points did not allow us to explore the region Aj/y > a, 
this is not expected to be physically interesting since the critical temperature and 
the fermion mass function is always orders of magnitude smaller than a. For the 
same reason, we tend to trust more our results with lower Ayy/a ratios, since 
then the lattice spacing becomes smaller and the integrations more accurate. At 
momenta of order Ajjv = 0.1a the mass function is already negligibly small, and 
thus our results are not influenced by such a UV cut-off. We were not able to 
consider even smaller A[/y/a as in || due to numerical problems caused by the 
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Fermion flavours =>■ 


N f = l 


N f = 2 


TV, = 3 


Af/y/a JJ- 


m 




r 


m 


t c 


r 






r 


0.1 


21 


4.5 


9.4 


3.3 


0.6 


11.0 


0.65 


0.12 


10.6 


0.5 


32 


6.2 


10.4 


4.0 


0.72 


11.1 


1.3 


0.26 


10.0 


1 


32 


6.1 


10.5 


4.4 


0.75 


11.7 


1.4 


0.27 


10.4 



Table 1: The quantities to = 10 3 x M(0,0)/a at T = 0, t c = 10 3 x k B T c /a and 
r = 2mo/t c for different values of flavours and ultra-violet cut-offs A^y. The vertex 
with Ta{j>, q) = 1 + A(q) is used. 



Fermion flavours =>■ 
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TO 
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0.1 


22 


4.6 


9.6 


3.4 


0.67 


10.1 


1.1 


0.21 


10.5 


0.5 


29 


5.7 


10.2 


4.7 


0.85 


11.1 


1.4 


0.27 


10.4 


1 


31 


5.8 


10.7 


5.0 


0.9 


11.1 


1.5 


0.28 


10.7 



Table 2: The same quantities as Table 1 using the symmetric vertex. 
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Figure 3: The fermion mass function at zero momentum and zero temperature, 
scaled by a and on a logarithmic scale, with respect to Nf for a ratio Auy/a = 0.1. 
Stars and squares correspond to the vertex with 1^4 = 1 + A(p) and the symmetric 
vertex respectively. We fit our results for both vertices with the curve M(0, 0)/a = 
e _1 " 5JV //13.1. Values of Nf larger than 3.35 are not considered, because then the 
self-energy falls below the IR-cut-off. 
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complexity of the algorithm. 

We note that in general the value of r remains roughly stable for various 
choices for Nf and A-uv- It is in fact more stable than the r ratio computed in ||, 
which could be an indication that the wave- function renormalization A(po, \p\) plays 
a -at least numerically- stabilizing role in the system under study. In this respect, 
the results presented here could be more trustworthy. Furthermore, comparing the 
Tables 1 and 2 one can see that the r-ratio is not considerably affected by the 
particular choice of our vertices. Moreover, we note that its value is concentrated 
approximately around r ~ 10.7, which is comparable to the values obtained in || 
which neglected A(jjq, \p\). The ratio r found is also comparable to typical values 
obtained in Ref. fl2fl , which includes wave- function renormalization effects but uses 
several approximations which we were able to by-pass in this study. We confirm 
therefore that adding these effects does not influence the behaviour of the theory in 
a significant way. We have to note moreover that our r values are somewhat larger 



than the value r w 8 measured for some high-temperature superconductors [IS]. 
However, as already discussed, we could be overestimating this ratio because of a 
possibly poor convergence of the algorithm for temperatures close to the critical 
one. 

3.3 Behaviour of the theory with the number of fermion flavours 

At zero-temperature, this theory is known to exhibit also an interesting behaviour 
with the number of fermions Nf. In Fig. 3 we plot the zero- momentum and zero- 
temperature fermion mass function with respect to Nf, fitted with the same ex- 
ponential curve for both vertices, since their behaviour is similar. There are two 



14 



reasons why we do not try to fit them with a function non- analytic in Nf, as was 
done in ||. One reason is that the non-analytic form used before and predicted by 
some studies is expected to be valid only for a number of fermion flavours close to 
a possibly critical value N c , so efforts to use the same functional form for Nf away 
from this value could give a misleading value for N c . The other is that our present 
results for the ratio r seem to be quite stable with Nf, and we would like to be able 
to estimate in a simple way some sort of "mean" r-ratio by combining these data 
with the ones showing the T c dependence on Nf, as will be described below. 

At Nf ~ 3.35, the mass function is still roughly four times larger than the cut- 
off. When Nf > 3.35, our algorithm does not converge and the mass function tends 
fast below the IR cut-off. Note that this value did not vary substantially for different 
choices of Am or the vertex. This behaviour could indicate that Nf 3.35 is some 
critical point beyond which dynamical mass generation is impossible. However, 
numerical limitations did not allow us to explore the region Am < 10 _4 Af/y. In 
other words, if there were a solution for the mass function smaller than roughly 
10~ 4 A[/y, our algorithm would not be able to find it because it would fall below 
the IR cut-off. The finite size of our lattice therefore does not allow us to draw firm 
conclusions on the matter, and the fitting curve should not be extrapolated beyond 
this limiting Nf value. 

The value of Nf where this behaviour comes in effect is remarkably close to 
the one quoted in the numerical study of ||, and it is also close to our previous 
result jaj. These values are however roughly 25% smaller than the one quoted in 



the theoretical studies of Refs. [20|, 21] which take wave-function renormalization 



effects into account. A similar study @, which includes a calculation of the fermion- 
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field anomalous dimension to second-order in 1/Nf, predicts a critical value Nf ~ 
3.28, which is only slightly larger than the theoretical prediction neglecting wave- 
function renormalization which gives Nf = f| w 3.24 0, and quite close to the 
value where our algorithm looses convergence. 

In Fig. 4 we plot the phase diagram of the theory with respect to Nf and 
JcbT. It separates two regions of the parameter space which either allow or do not 
allow dynamical mass generation. The approximate stability of r with increasing 
Nf lead us to rescale the data points of this plot and fit these simultaneously with 
the data points of Fig. 3, with the same exponential curve e~ 1 ' 5N f /13.1 for both 
vertices, since they exhibit a similar behaviour. The slope of this exponential fits 
also relatively well the data on the two plots separately, even though the slope of 
the exponential of the mass function separately is closer to -2. It is however lower 



in absolute value than the one reported in [15] . We note also that the fall of the 
critical temperature with the number of fermion flavours is slower than in ||, which 
was found to be -2.5. The fast fall in that reference was the cause of the increase of 
the r-ratio with Nf , and we interpret it as poor convergence of that algorithm with 
large Nf values, because then the self-energy becomes smaller and more vulnerable 
to instabilities coming from numerical errors. The (double of the) rescaling factor 
was found to be f = 2M(0, 0)/fcgT ~ 10.7 and provides some kind of "mean" r- 
ratio which is independent of Nf. Its approximate value could actually be guessed 
by the data in Tables 1 and 2. 

The choice of an exponential fitting curve was only made to describe "phe- 
nomenologically" the general tendency of the data and to provide a measure for a 
r-ratio independently of Nf, and is reminiscent of the results in Ref. || but with a 
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Figure 4: The phase diagram of the theory with respect to the temperature scaled by 
a and in logarithmic scale and the number of fermion flavours, for a ratio Ajjy /a = 
0.1. Stars and squares correspond to the vertex with Tjy = 1 + A(p) and the 
symmetric vertex respectively. We fit the data points with the curve ksT/a = 
e -i.5N f / 70 _ Thig curye shoulc j not be ex trapolated for N f > 3.35. 
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steeper slope. However, there are also studies that predict a non-analytic behaviour 
of £ for Nf near its critical value Q. One could speculate that such a behaviour 
comes into effect also for our larger Nf data, since the slope of a curve connecting 
the large-iVy points of both Fig. 3 and 4 is larger than the one of the general fit. 
Moreover, as already noted, near these Nf values the mass function is still roughly a 
factor of four larger than the IR cut-off, so the value of Nf quoted cannot be easily 
argued to be just an artifact of the finite lattice size. Lack of convergence of the 
algorithm and fall of the mass function below the IR cut-off however do not allow 
us to test the precise behaviour of the theory near Nf ~ 3.35, and it is also clear 
that the fitting curves should not be extrapolated for Nf larger than this value. 

We remind the reader that the results quoted above were obtained with a 
photon propagator in the Landau gauge, which is convenient and widely used for 
this type of calculation. In addition, it possibly enables the wave-function renor- 
malization to play a stabilizing role in the numerical algorithm and provide more 
accurate results. Since the final answers should be gauge-invariant however, one 
should in principle check whether other gauges, with corresponding vertices con- 
sistent with the Ward-Takahashi identities, give the same answer. In particular, 
calculations in the non-local gauge |22| eliminate the wave-function renormalization 
dependence of the results, since A{p) = 0, which means that the ratio r or the 
possibly critical number of fermion flavours should not be influenced by A{p). In 
such a case, however, a different vertex choice than ours should be made, so that 
it does not reduce to the bare one when A(p) = 0. In our study, therefore, it is 
the taking into account of the wave- function renormalization in conjunction with a 
non-bare vertex that could potentially lead to results different from || . 
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4 CONCLUSIONS 



To conclude, we were able to solve a system of two integral equations for the fermion 
mass function and wave-function renormalization for a finite-temperature version 
of three-dimensional QED, by applying a numerical relaxation technique. The 
results presented are of course specific to the dimensionality of the system under 
study and to the abelian character of the theory. Both functions as well as the 
photon self-energies are taken to be energy- and momentum-dependent, use of two 
different non-bare vertices is made, and this leads to a more accurate and detailed 
study of the criticality of the theory. One main result is a r-ratio of about 10.7, 
which is close to previous numerical studies, confirming that including wave- function 
renormalization, in conjunction with two different non-bare fermion-photon vertices, 
does not affect the theory in a significant way. The other important result is the 
existence of a possibly critical fermion flavour number of roughly 3.35, which is 
also consistent with some theoretical expectations and other numerical calculations. 
But the use of an IR cut-off does not allow the precise exploration of the possible 
criticality at this stage. 

It should be stressed here that such an agreement could not be taken for 
granted a priori, in view of the seriousness of truncations and approximations used 
in the past. It is the first time that such a study drops most of the approximations 
applied so far and allows a reliable description of the behaviour of the theory. We 
estimate the numerical uncertainty for the values quoted, which comes mainly from 
the convergence criteria imposed, at about ±10%. The specific truncations of the 
Schwinger-Dyson hierarchy chosen in this work could of course easily introduce an 
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additional theoretical error of comparable magnitude. We also note that a larger 
lattice than the one we used could be expected to improve the convergence of the 
algorithm and help us understand the behaviour of the fermion mass function at 
regions even closer to the critical temperature and flavour number. Furthermore, it 
remains to be seen if the inclusion of the imaginary parts of the self-energies could 
influence these results considerably. 
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